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LARGE  AMPLITUDE  MULTIMODE  RESPONSE  OF  CLAMPED 
RECTANGULAR  PANELS  TO  ACOUSTIC  EXCITATION 


By 

Chub  Mei 

INTRODUCTION 


Acoustically  induced  fatigue  failures  in  structural  components  have 
resulted  in  unacceptable  maintenance  and  inspection  burdens  associated  with 
aircraft  and  missile  operation.  In  some  cases,  sonic  fatigue  failures  have 
resulted  in  major  structural  redesigns  and  aircraft  modifications.  Thus, 
accurate  prediction  methods  are  needed  to  determine  the  fatigue  life  of 
structures. 

Many  analytical  and  experimental  programs  to  develop  sonic  fatigue 
design  criteria,  however,  have  repeatedly  shown  a  poor  comparison  between 
measured  and  calculated  maximum  RMS  stress/strain  (refs.  1,  2).  Deviations 
in  excess  of  100  percent  are  not  uncommon.  Large  deflection  nonlinearity 
has  been  identified  as  a  major  factor  for  the  enormous  discrepancy  between 
test  data  and  computed  results  (ref.  3).  A  test  program  was  conducted  re¬ 
cently  to  check  the  analytical  effort  for  the  large-ampl itude,  single-mode 
response  reported  in  reference  3.  The  acoustic  response  tests  were  perform¬ 
ed  in  the  Wideband  Acoustic  Facility  at  Wr ight-Patterson  Air  Force  Base.  A 
comparison  of  the  results  from  two  panels  is  shown  in  figure  1.  The  predic¬ 
tion  of  random  responses  is  much  improved  with  the  single-mode  computational 
method,  especially  at  high  excitation  levels.  Test  results  (fig.  2)  also 
showed  that  there  are  more  than  one  mode  responding.  Multiple  modes  were 
also  observed  by  White  in  experimental  studies  on  aluminum  and  carbon  fiber- 
reinforced  plastics  (CFRP)  plates  under  acoustic  loadings  (ref.  4).  White 
also  showed  that  the  fundamental  mode  responded  significantly  and  contribut¬ 
ed  more  than  80  percent  of  the  total  mean-square  strain  response;  higher 
modes,  up  to  third  or  fourth  modes,  account  for  95%  or  more  of  the  total 
mean-square  strain  response.  In  order  to  have  an  accurate  prediction  of  the 
random  response  of  a  structure,  multiple  modes  should  be  used  in  the 
formulation. 


RESULTS  COMPARISON 


Figure  1.  Comparison  of  analytical  and  experimental  mean— square  stresses  of 
clamped,  square,  aluminum  panels. 
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AMPLITUDE  MICRO-STRAIN 


NOMENCLATURE 


a,b 

C 

D 

E 

fm(x) >Sn(y) 

F 

h 

K 

L 

M 

P 

q 

r 

Sp(u) ) 
t 

w 

W 

x.y 

6 

Z 

0) 

8 

Subscripts 

EL 

L 


panel  length  and  width 
generalized  damping 
flexural  rigidity 
Young* s  modulus 

displacement  functions,  eq.  (21) 

airy  stress  function 

panel  thickness 

generalized  stiffness 

mathematical  operator,  eq,  (1) 

generalized  mass 

pressure 

normal  ccoordinate 
length-to-width  ratio,  a/b 
cross-spectral  density  of  p(t) 
t  ime 

lateral  deflection 
generalized  displacement 
coordinates 

vector  function,  eq,  (16) 
damping  ratio,  c/cQ 
normal  mode 
linear  frequency 

equivalent  linear  or  nonlinear  frequency 

equivalent  linear 
linear 


MATHEMATICAL  FORMULATION  AND  SOLUTION  PROCEDURE 


The  governing  equations  of  a  rectangular,  isotropic  plate  undergoing 
large-deflection  motions,  neglecting  the  effects  of  both  inplane  and 
rotatory  inertia  forces,  are  (refs.  5,  6) 

L(w,F)  =  DV^w  +  Phw»tt  +  gw»t 


-  h(F ,  w, 

yy 


XX 


F,  w, 
xx  yy 


-  2  F,  w 

xy 


,  ) 
xy 


-  p(t)  =  0 


(1) 


V^F  *  ECw2, 

xy 


W, 


w 


XX 


•yy> 


(2) 


where  a  comma  denotes  the  partial  differentiation  with  respect  to  the 
corresponding  variable,  w  is  the  lateral  deflection,  F  is  the  stress 
function,  D  is  the  flexural  rigidity,  P  is  the  mass  density,  h  is  the 
plate  thickness,  p  is  the  pressure,  E  is  the  Young* s  modulus,  and  g  is 
the  viscous  damping. 

The  lateral  deflection  is  assumed  as 


w(x,y,t)  =  h  l  l  W  (t)  f  (x)  g  (y)  m,n  =  1,2,3,... 
mn  m  n 

m  n 


(3) 


where  the  functions  ^(x)  and  gn(y)  are  so  chosen  that  they  satisfy 
the  boundary  conditions.  By  solving  the  compatibility  equation,  equation 
(2),  the  stress  function  can  then  be  determined  as 


F  =  N  £-  +  N  ^ 
x 


f-  +  \  -y  +  Eh2  l  I  F..  N^x)  M.(y) 


i,j  »  0,1,2,. 


(4) 
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A  quasi-exact  solution  has  been  obtained  by  Paul  for  thermal  postbuckling  of 
a  clamped,  rectangular  plate.  The  expressions  for  the  coefficients  N  , 

N  ,  and  F  can  be  found  in  reference  7. 

y  ij 

Apply  the  Bubnov-Galerkin  method  to  the  equation  of  motion  in  deflec¬ 
tion,  equation  (1),  as 

//  L(w,F)  f  g  dxdy  =  0  r,s  =  1,2,3...  (5) 

r  s 

After  performing  the  integration  over  the  total  area  of  the  panel*  a  set  of 
nonlinear,  time-differential  equations  is  obtained  and  can  be  written  in 
matrix  form  as 

[M]{W>  +  [C]  {W>  +  [K]l  {W}  +  (3  (W)}  =  { p( t)}  (6) 

where  the  matrices  [M] ,  [C],  and  [K]l  are  the  generalized  mass,  damp¬ 
ing,  and  linear  stiffness  matrices,  respectively,  and  { 3 (W ) }  is  a  vector 
function,  cubic  in  the  generalized  displacements  {W} . 

An  equivalent  linear  set  of  equations  to  equation  (6)  may  be  defined  as 
(refs.  8-13): 

[M]  {W}  +  [C]  (W>  +  (  [K]  L  +  [K]el)  {W}  =  {p(t)>  (7a) 

or 

[M]  {W>  +  [C]  {W}  +  [K]  {W>  =  { pCt)}  (7b) 

[K]  can  be 

Ci  It 

(8) 

where  E[  ]  is  an  expected  value  operator. 


where  the  elements  of  the  equivalent  linear  stiffness  matrix 
obtained  from  the  expression 


"33 ,  (W) 

~T — 
ij  3W. 


i,j  1)2,3,... 
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To  determine  the  mean-square  generalized  displacements  F^  in  equation 
(7),  an  iterative  solution  procedure  is  introduced.  The  undamped  linear 
equation  of  equation  (7a)  is  solved  first.  This  requires  the  determination 
of  the  eigenvalues  and  eigenvectors  of  the  undamped  linear  equation 

0)2  [M]  {<J>K  =  [k]l  UK  (9) 

where  ai j  is  the  frequency  of  vibration  and  {(J)}j  is  the  corresponding 
normal  mode  shape  based  on  linear  theory. 

Apply  a  coordinate  transformation,  from  the  generalized  displacements 
to  the  normal  coordinates,  by 

{W}  =  [<H  {q}  n  <  m  (10) 

mx  1  .m*n  n*l 


in  which  each  column  of  [$]  is  a  modal  column  of  the  linear  system,  and 
{  q)  represents  the  normal  coordinates.  Substituting  equation  (10)  into  the 
damped  linear  equation  of  equation  (7a)  and  premultiplying  by  the  transpose 
of  [<f>  ] ,  it  becomes 

[m]  {q}  +  H  {q}  +  [k]l  {q>  =  (P(t)}  (11) 

where  =  t<)>]T  [M]  [<(>  ] 

[k]l  =  [$]T  [K]l  [*]  =  Qo2]  [m] 

[c]  =  [4»]T  [C]  [<(.]  =  2^]  M 

{P}'  -  [<MT  (p>  (12) 


The  mean-square  normal  coordinate  is  simply 


(14) 


irS  (10.)  ir{ <}>} .  [S  (o>.)]  {<()} 

_ 2 _ J _  _ _ 2 _ P_J _ ± 


4  r  «d? 
J  3  J 


4  M?  Z. 

J  J  J 


where  [S  ]  is  the  cross-spectral  density  matrix  of  the  excitation  { p( t ) } . 
P 

The  covariance  matrix  of  the  linear,  generalized  displacements  is 


w{*}.  [S  (\)1  4>k  T 

[WjW.J,  -  I  U>k - k  !  - . Wk 

1JL  k 


The  diagonal  terms  [W^W^ ] ^  are  the  mean-square,  linear,  generalized  dis¬ 
placement  W2.  This  initial  estimate  of  W2  can  now  be  used  to  compute 
the  equivalent  linear  stiffness  matrix  [K]^  through  equation  (8).  Then 
equation  (7)  is  again  transformed  to  the  normal  coordinates  and  has  the  form 


M  +  [c]  {q}  +  |>]  {q}  =  {P(t)J 


where  [k]  =  [4>]T  ([k]l  +  [K]^)  [+]  =  [a2]  Cm] 


The  jth  row  of  equation  (17)  is 


q  +  2Cj  q^  +  a2  q  -  ^ 


and  the  displacement  covariance  matrix  is  given  by 


_  »(♦>£  [S  (SL )]  4>k  T 

[W  w  ]  -  l  {*}  -  p  - -  4>k 

J  k  4  “k  «k  \ 
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Convergence  is  considered  achieved  whenever  the  difference  of  the  RMS 
generalized  displacements  satisfies  the  requirement 


(RMS  W,).  "  (RMS  W.).  . 

j  lter _ J  lter~1  <  10-3,  for  all  j 

(RMS  W.). 

j  iter 


Once  the  RMS  displacements  are  determined,  the  RMS  deflection  of  the  panel 
and  the  maximum  RMS  strain  can  be  determined  from  equation  (3)  and  the 
strain-displacement  relations,  respectively. 
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DEVELOIMENT  OF  GENERALIZED  MATRICES 
AND  COMPUTER  PROGRAMS 


The  deflection  of  the  panel  is  represented  by 


w(x,y,t)  =  h  l  I  Wmn(t) 
in  n 


cos 


(m  -  1)ttx  (m  +  l)irx 

-  -  cos  - 


(n  -  l)iry  .  _  _  (n  +  l)xy 
cos  - — ^ ~  cos  - £ - - 


(21) 


The  expression  w  satisfies  the  boundary  condition  for  clamped  edges: 

w  =  w,  =0  on  x  =  0  and  a 
x 

w  =  w,  =0  on  y  a  0  and  b  (22) 

y 

The  stress  function  can  be  expressed  in  terms  of  the  generalized  dis¬ 
placement  as 


F  -  \  £  +  \  T  +  l  *  Fij  cos  TT  cos 


211 

b 


(23) 


where  the  coefficient  F-jj  is  given  by  the  expression 


«  /i* 


(r- +  j2r) 


-  I  l  l  l  B.  .  .  -  W  W, 

L  l  l  l  ijmnkl  mn  kl 
z  m  n  k  1 


(24) 


in  which  B.  .  , ,  are  integers  and  r  =  a/b.  The  coefficients  N  and  N  , 

ljmnkl  x  y 

and  the  integers  are  given  explicitly  in  reference  7.  The  partic¬ 

ular  generalized  displacements  that  are  chosen  to  be  nonzero  in  the  conver¬ 
gence  studies  are  shown  in  table  1. 
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Table  1. 


Generalized  displacements  for  convergence  studies. 


Number  of  terms 


Generalized 

Displacements  4  6  _10  _15_ 


wn 

Wl3 

W31 

W33 

W15 
W51 
W3  5 
W5  3 
W17 
W71 
W55 
W37 
W7  3 
W19 


W( 


91 


XXX  X 

XXX 

XXX 

XXX 
X  X 

X  X 

X 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 
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Utilizing  the  expressions  for  w  and  F,  equations  (21)  and  (23), 
respectively,  and  performing  the  integration  of  equation  (5),  the  integral 
associated  with  the  inertial  force  term  in  equation  (1)  has  been  derived  as 


fh  J*  Phw,tJ;  fr  gs  dxdy  -  (wr_2,s-2  -  2  Wr_2>s  -  2  WrjS_2 

o  o  \ 


+  4  W  -  2  W 


-  2  W 


r,s  r,s+2  r+2,s 


+  w  +  W  +  W 

r-2,s+2  r+2, s-2  r+2,s+2i 


(25) 


The  generalized  mass  matrix  [M]  in  equation  (6)  using  15  terms  in  the 
deflection  function  is  given  by: 


W11  W13  W31  W33  W15  W51  W35  W53  W17  W71  W55  W37  W73  W19  W91 


(26) 
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A  subroutine  program  MASS,  which  generates  the  mass  matrix,  has  been  coded 
and  verified.  A  listing  of  the  MASS  subroutine  is  given  in  the  Appendix. 

Similarly,  the  integrals  associated  with  the  linear  stiffness  terms  in 
equation  (1)  yield 


fb  ra  _34w  _  DhiT4ab 

/  J  D  —  fr  gs  dxdy  =  7“ 

o  o  3x4  4a4 

*  {  [  (r  -  l)4  +  (r  +  1)-J  [  (Ci  +  »Ht>,  -  Wt_s.2  -  WIjl+2] 

*  <r  -  1)"  [«r_2',_2  ♦  »r.2>8*2  -  (Cl  *  UVlJl 

+  (t  *  1)“  IWr.2,s*2  *  Wr+2,s-2  ‘  (C1  *  U  "r*2,  J) 

rb  ra  34w  _  ,  _  Dh7r4ab 

J  J  D  f  gs  dxdy  =  — 

o  o  3  y4  4b4 

*  {  [  (s  -  l)1*  +  (s  +  l)4]  [  (C2  +  l)wr>s  -  Wr_2>s  -  Wr+2)S] 

+  (•  -  D4  [wr-2,s-2  +  wr+2, s-2  “  + 

+  <s  +  D4  CWr+2,s+2  +  Wr-2,s+2  ’  <C2  +  Wr,s+2]} 


fa  2D  _ !L-^ _  f  g  dxdy 

00  3k2  ay2  r  S 

*  {  (r  —  l)2  (s  +  l)2  [w 

r,£ 

+  (r  -  l)2  (s  +  l)2  [w 

r,s 

+  (r  +  l)2  (s  -  l)2  [w 

r,  s 

+  (r  -  l)2  (s  -  l)2  [w 


Dhir4  ab 


-  W  -  W  +  W  J 

r, s+2  r+2, s  r+2,s+2 

-  w  —  w  +  w 

r,s+2  r-2,s  r-2,s+2 

'  "t,.-2  -  Wr.2,s  +  V2,.-2] 

’  “t,s-2  •  Wr-2,  s  *  V2.S-211 
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The  generalized  linear  stiffness  matrix  [K]l  equation  (6)  using  15 
terms  in  the  deflection  function  has  been  derived.  It  can  be  expressed  as 
the  sum  of  the  three  submatrices  as 

[K]  =  Dh--—  ( —  [K] !  +  —  [K]2  +  [K]  3^  (30) 

L  4  \  a4  b4  a2b2  ) 

The  nonzero  elements  of  the  three  linear  stiffness  submatrices  are  given. 
Since  the  stiffness  matrix  is  also  symmetric,  only  the  lower  left-hand  side 
elements  are  given.  They  are 


W11 

W13 

W31 

W3  3 

W15 

24(Ci+l) 

-24 

24(C1+1) 

? 

-2^  (Ci  +1) 

24 

(24+44)(Cl+l) 

24 

-24 (C1+i) 

-(24 +44) 

(24+44)(Ci+1) 

0 

"24 

0 

24 

24  (  C  l +1 ) 

K]l  = 

0 

0 

-44(C1+1) 

44 

0 

0 

24 

0 

-(24+44  ) 

-24(Ci+1) 

0 

0 

44 

-^(C^l) 

0 

0 

0 

0 

0 

-24 

0 

0 

0 

0 

0 

W51 

W3  5 

W53 

W17 

W71 

^+6I+)(C1+1) 

0 

(24+44)(Cl+l) 

symmetric 

-(44  +64 ) 

44 

(44+64)(C1+l) 

0 

24 

0 

24 (c 1+1) 

+1) 

0 

64 

0  (64+84  )(C1+1)J 

(31) 
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W51  W35  W53  W17  W71 


24  (C2+l) 

0 

(44  +64)(C2+l) 

symmetric 

-  2ft  (  C2  +1 ) 

44 

(24 +44 ) (C2+l) 

0 

64 

0 

(64+84)(Ci+1) 

0 

24 

0  24(C!+1) 

(32) 
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W31 


[k]3 


W11  W13 

(2*  2)2 

-(2*  2)2  (2 •  4) 2  +( 2 * 2)2 

"(2*  2)2  (2*  2)2 
(2*  2)2  — ( 2 •  4-) 2  — ( 2 •  2) 2 
-(2*  4)2 
0 

(2*  4)2 
0 
0 
0 


0 

0 

0 

0 

0 

0 


W33 


w 


15 


(4*  2)2  +(2*  2)2 

-(4*2)2-(2*2)2 

0 

-(4*2)2 

0 

(4*2)2 

0 

0 


symmetric 

(4*4)2 

+2(4*2)2+(2*2)2 

(2*4)2 

(4*2)2 

-(4*  4)2  -(2*  4)2 
-(4*4)2-(4*2)2 
0 
0 


(2*  6)2 +(2*4)2 
0 

-(2*6)2-(2*4)2 

0 

‘(2*6)2 

0 


W, 


51 


W- 


3  5 


W, 


53 


W 


17 


W 


71 


(6*  2)2+(4*  2)2 

(4*  6)2+(2*  6)2 
0  +(4*  4)2+(2*4)2 

(6* 4)2+(4*4)2 

-(6*  2)2-(4*  2)2  (4*4)2  +(6*  2)2+(4*  2) 2 

0  (2*6)2  0  (2*  8)2+(2*  6)2 
-(6*2)2  0  (6*2)2  0  (8*  2)2+(6*  2)2 

(33) 


The  nonzero  elements  of  the  linear  stiffness  matrix,  kL(i,j)  for  i,j  >  11, 
are 


kjO.1,4)  =  kx  (12, 11)  =  44 

kx (12, 5)  -  kx (14,12)  =  24 

kx  (13, 6)  =  kx  (13, 11)  =  64 

kx  (11,7)  =  -44(C1  +  1) 

k,  (12,7)  =  ~(24  +  44 ) 
k^ (11,8)  =  -(44  +  64 ) 


(34) 

(cont 'd) 
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kx  (13,8)  =  -64  (Cx  +  1) 

k: (12, 9)  =  -24  (Cx  +  1) 

kx  (14,9)  =  -24 

kx  (13,10)  =  -(64  +  84) 

kx (15,10)  =  -84  (Cx  +  1) 

kL  (11,11)  =  (44  +  64)(C1  +  1) 

kx  (12,12)  =  (24  +  44)(C1  +  1) 

kx (13,13)  =  (64  +  84)(C1  +  1) 

k: (15 , 13)  =  84 

kx  (14,14)  =  24(Cx  +  1) 

kx  (15 , 15)  =  (84  +  lCh)(Cl  +  1) 

1*2(11,4)  =  1^(13,11)  =  44 
1*2(12,5)  =  1*2(12,11)  =  64 
1*2(13,6)  =  1*2(15,13)  -  2f* 

1*2(11,7)  =  ~(44  +  64 ) 

1*2(12,7)  =  -64(C2  +  1) 

1*2  (11,8)  =  -44  (C2  +  1) 

1*2(13,8)  =  -(24  +  44) 

1*2(12,9)  =  -(64  +  S4) 

1*2(14,9)  =  -8tt(C2  +  1) 

1*2(13,10)  -  -24(C2  +  1) 

1*2(15,10)  =  -24 

1*2  (11,11)  =  (4  +  6)(C2  +  1) 

1*2  (12,12)  «  (64  +  84)(C2  +  1) 

1*2(14,12)  =  84 

1*2(13,13)  =  (24  +  44)(C2  +  1) 

1*2(14,14)  =  (84  +  10f»)(C2  +  1) 

1*2(15,15)  =  24(C2  +  1) 

k3  (11,4)  =  (4‘4)2 

k3  (12,5)  =  k3  (13,6)  =  (2‘6)2 

k3 ( 1 1 , 7 )  =  k3 (11,8)  =  -(4*  4)2  -  (4*  6)2 

k3 (12,7)  =  k3 (13,8)  =  -(2.6)2  -  (4.6)2 

1*3(12,9)  =  k3  (13,10)  =  -(2-  6)2  -  (2*8)2 

k3 (14, 9)  =  k3 (15, 10)  =  -(2.8)2 

k3  (11,11)  =  (6*6)2  +  2(6*4)2  +  (4»4)2  (34) 

k3  (12,11)  =k3(13,ll)  =  (4*6)2  (cont'd) 
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k3  (12, 12)  =  k3 (13,13)  =  (4*8)2  +  (2*8)2  +  (4*6)2  +  (2* 6) 2  (34) 

k3 (14,14)  =  k3 (15,15)  =  (2.8)2  +  (2.10)2  (concl  d) 


where 


(  for  s  -  1 
(forsM 

f for  r  =  1 
I  for  r  *  1 


(35) 


A  subroutine  program  LSTF  which  generates  the  linear  stiffness  matrix  has 
been  coded  and  verified.  A  listing  of  the  LSTF  subroutine  is  presented  in 
the  Appendix. 

Derivation  of  the  generalized  equivalent  linear  stiffness  matrix  [K] 

EL 

in  equation  (7a)  has  been  initiated.  It  is  in  good  progress.  Continuing 
research  effort  will  be  devoted  to  the  following  tasks: 

(1)  Completion  of  the  derivation  of  equivalent  linear  stiffness; 

(2)  Application  of  eigen  solution  and  coordinate  transformation; 

(3)  Determination  of  mean-square  linear  generalized  displacement; 

(4)  Implementation  of  the  iterative  process; 

(5)  Derivation  of  strains  computation; 

(6)  Coding,  debugging,  and  verifying  the  complete  computer  program; 

(7)  Convergence  studies;  and 

(8)  Generation  of  design  charts. 7 

o 
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APPENDIX 


LISTINGS  OF  THE  MASS  AND  LSTF  SUBROUTINES 


f 

ji  '  mass 

A 

LN46O0.FOR  FORTRAN  V.5(515)  /KI  29-JUN-81 

14138  PA(|E  1 

/H  aeesn 

- S«TgR (JTJTTN E'  MASS  (  AO#  FL^H,  RHO,  5M,  NTERM  j 

■2100602 

DIMENSION  SM(NTERM,NTERM) 

j  3  j 

C 

Ui  09604 

c 

'  THIS  SUBROUTINE  GENERATES  THE  SYSTEM  MASS  MATRIX  OF 

THE  PANEL-  USTNTG 

j  s ! 06005 

c 

CNTEFtM)  TERMS  IN  THE  DEFLECTION  FUNCTION 

■  - 

jei  00086 

c 

ALaPANEL  LENGTH - 

I  / 1  8***3  £2  7  . 

-~C~ 

•  SL=¥S¥£L'  width 

{■j  109098  < 

c . 

HsPANET,  THICKNESS 

..  '  .  :-i  /.  '  "■ 

j  s I  000^9 

c 

RHOaMASS  DENSITY 

; -si  k}$Y}  1  :>* 

~~ (T 

Sh  ( frrfchM-;  NTERM)  sSTsTEW*^.  "CEllEH7fCTZE!T"ff A £5~ MATR fX 

;i k>0k)  1  1 

C 

WTErmsI#  4 1  6#  10#  Or  15 

j  — •  'j&tt  1  2 

C 

i  '*3  »£0v?  1  3 

CTOTae  .  25  #ft'HO*K*H*lE*RE 

*  }i4f  J4 

. c 

INITIALIZED  THE  MASS  MATRIX 

K08wl5 

00  10  I®1  #NTErM 

v'  :  -  ..  ‘  ‘  : 

it  6:  000  1  O 

DO  10  Ja  1  #  n’TEPm 

h'j  7 

SM(  I#  J)=(;  .3 

|l3:  00018 

10  CONTINUE 

|  I3.  00019— 

- S^f  1#  H  =  4  .0*COFF 

•  - 

IF  (NTEHM  *EQ,  t)  GO  TO  26 

•  4/  -  4:;.  v  ’  • 

*21*  00021  .; 

WW 

1  v:. 

M 00622 

SM(1#2)=-2.0*COeF 

izjj  yd*?  2  3 

SM(l,3)=-2.0#C0EF 

jaMj  0002  4 

SM ( 1 »4 ) aCOEF 

S^(2f2)34t0*COEF 

•as!  00026 

.  SM(2#3)sCOEF 

..  ..  . 

M 00027 

SM  (  2  »  4  )  -*»2  .*?*C  0  EF  . 

!  28100028 

j2»  @0629 

1  5M( 3  #4)^*2 •0*COEF 

'so!  03030 

SM ( 4  »  4 ) a4 ,0*COEF 

'31  (  *>«>¥/ M 

IF  C w TERM  m bU *  4}  60  TO  20 

i32:  000  3  2  > 

c 

.•  •:/>;-  : -/•  .-A  * 

■32:  ^0031  , 

-  '  sH(2rS)*-2.0*COFF 

!•  ' 

j24i  kj'^k)  3  4 

- SM(  3  ,<)  rrs-Ti^CCrrF  ■  “  '  ““ 

;js.  O003S 

SM(4#5)=C0EF 

be! 

S*(4,6)=C0EF 

;37  IT 

~  5w(5#5  7-4^(y»CuEF  . ...  ...  .  •  \  ■■  -  . ,.  T-  ' 

;  smc6#s)*4.^*coef 

•v;  ,  ••  •;.-  ;• 

]39!  000.J#>' 

IF  (NTERM  ,E0.  6)  GO  TO  2* 

•;..:yv;  .>:y"  "V ' 

jaof  00049  . 

(T 

|4t !  000  4 1 

SM(2#7)aC0EF 

U!  08042 

SM(3#8)aCOEF 

zzrmoTTr 

- - 5«t4TTT=»2v0*COtF— ~ 

i,{4  90044 

sM(4#R)s-2.0*COEF 

•.«»;  60045;.;. 

".  S  M  (  5  *  7  )  a  •  2  »  0  *C  o  E  F 

;  ’  '•  :  ••  ’  '  '•  •' '  •  . 

irto!  6064B 

5^3777^2 . 0  *COTF  “  • 

1^00047 

SM(6»8)a-2.0*COEF 

j48i  00048 

SM(6#10)a-2.0*COEF 

iJ  00048 

SM  (7  #7  )  *4 ,0#COEF - " — "  1  .  -  " 

..V  -1  •• 

-'so!  00050  J 

SM(7,8)aCOEF  . 

Jsii  0805 1  < 

SHC7>9)aC0EF 

V  V 4;.:  ...  '  '  ‘  -  • 

]53i  00053 

SM(8#10)aCOEF 

1541  00054 

1  _ 

SM(9»9)s4.0*COEF 

;§3i  MJ.35."..'. 

$M  (  10,  1  !7)tf4‘.0*COEF  '  ■ 

jsei 

IF  (NTERM  ,EQ.,  10J  GO  TO  20 

... 

1  11  . . 1  11  . . . . . . . . . . 1"" . . . .  " .  1 . 
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HASS 


LN4PQ0.FOP 


FORTRAN  V.5(515)  /Kl 


29-iJUN-8i 


14*38  PAGE  1 


!a;  *s»058  SM(4, 1 l )sCOEF 

Js!  ; :.U  5^C5#  12)*C0£F  . 


!  4;  MVO*  S’  lb#  1  j;3LUC.f 

!  si  00061  SMC7, 11)s-2.0#C6EK 

1^00262  5M(7,12)s-2.0*COEF 

iThmrwi - - sin'T7ii'Ts-2.0*coEr  ” 

i  a !  00064  SM<3#13) *«2 *0*COEF 

Uiw^wfeS  SM(9,12Js-2.0*CoEP 

'—-nr^Fb - — — TT7#“r4Ts^?7r*cr>  ff~ 

iui  20067  'SM(  ia»  13)*-2.w*COEF 

00*6 »  SM  ( tv- #1*$)  «»2  .  n*r  OFF 

< r-si  wpwsr — — ■ — -  s^fii  #  i  V vs* . s »coef'  ' 
■i4  amlZr  ,  sM(u»  12)*coe:f 

iiaflMTI-  SMU1#13)»C0EF 

. . . s>Tf2ii2 ) «TT0**ci rzr^ 

-•it!  00073  SM(12#  14)*COEF 

!,3:»0074  SH  ( 1  3  #  13 )  s4  ,id*COEF 

TsmJWTT - — 1371  STsCOET  -  - — - 

iaoi  00076  :  5M  C  l 4 #14)  >4  ,0*COEF 

’sii  00077  ;  S  SM(15#15)a4,0*COEF 


y 00079  00  3«  Jsi,NTERM 

y 00080  DO  30  IaJ,«TEBH 


•IS’  VOV  ~  *  J»n  l  *  #  U  /  l  u  #  J,  J 

!m!  0008  2  *' :  3D  CONTINUE  .* 

;«>  0008  3  '..  .  •  ,  - :  RETURN 


Issi  SCALARS  AMO  ARRAY'S  C  «*«  MO  EXPLICIT  DEFINITION  .  »%»  NOT  REFERENCED  J 


30;  *<J  S 

;38,  »RHO  ;;.:y  13 

i«oi  •  l«5i}i30  50 

Li! 

Us!  temporaries 


»COEF — "  7~ 
v.S0003  7 
*AL  .  1,4 


•*H - - -  ~"3 

.50002  10 
*10002  15 


.SP001  11 

*1,0001  t&n 


DM - 5— 

.seem  12 

*1  ■  17 


•  / — - -  ■  -  ■  ■  ;  ~~ 
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c 

Li 

mm  1  -  .  SUBROUTINE  IjSTF (  AL#  Bl»#  H#  0  *  PI  # SK# NTERM J 

mm2  DIMENSION  SK(NTERM,«TERM) 

viniivM'  r  .  .  . 

I; 

Li 

mm*  C  THIS  subroutine  generates  the  system  linear  stiffness  matrix 
mm*  C  USING  (NTERM)  TERMS  in  the  deflection  function  for  the  panel 

03036  C  ALsPANEL.  LENGTH 

1: 

I  9 

mrnl  C  5)L»PA#»Et  “IDTH 

mm 8  c  hspukel  thickness 

<*00«9  e  OaE'*R**7/f  ft  7.  a*H-V*#7llsBENDrRG  RIGIDITY  . 

[To 

tt 

12 

30010  c  v=roisscr.,#s  ratio 

30011  ,  C  3K(NTERM,NTERM)rSYSTEM  OR  GENERALIZED  STIFFNESS  MATRIX 

l  2  C  UTERf-'.al  •  4,  6,  U .  OSL  IS 

71 

14 

13 

0301 3  C 

0S014  PI4aPI*PI*PI*Pl 

r  o  p  p=<i  T  *  n  *  w*pt  4  *  h  r,  *  ht.  . 

16 

17 

18 

00016  A4sl.0/(AL*AL*AL*AL) 

30017  B4sl ,0/(BL*BL*BL#BL) 

30018  AB2s4.0/(AL*AL*BL*BL) 

19 

20 

21 

00019  C  ClsCl  FOR  S<>l,Cl=CCl  FOR  S*1 

30023  C  C2sC2  FOR  R<>t,C2»CC2  FOR  R»1 

0002 1  Ctsl-tf  .  s; ^ 

22 

23 

24 

00022  C2*l .0  •  > 

03023  CGls2.0 

00024  CC2a2T0 

23 

26 

27 

00025  C  INITIALIZED  THE  STIFFNESS  MATRIX 

00026  DO  10.  Iaf ,NTERM 

00027  no  u  Jal  f.«>rFRM  •.  . 

28 

28 

30 

00028  SK(I,J)=0.0 

30029  10  CONTINUE 

00030  SKM  »11sf  fCCt+1  t01*A4+(CC^  +  l  .0)*B4+AB2)*COEF 

iT 

32 

33 

0001 t  IF(MTERM  .EG.  1)  GO  TO  20 

00032  C 

000  1 3  SNf  1 , 21  =  -f  fCC2*l  .:O*B4  +  A4+AB21*C0EF  . 1 . S; J . 

34 

33 

1 

36 

03034  SK(l,3)s-((CCl+l ,0)*A4+B4+AB2)*COEF 

O0035  SK(l,4)a(A4+B4+AB2)»COEF 

013036  SK(7,2)  =  (  tT.P^rCC^  +  l  .i)I#B4  +  (Cl  +  l  .0  I*  A4+5  -fl»  AB2l*C0EF 

37 

38 

39 

000  37  3Kr2,3}*DA4+H4+AB2)*C0FF 

00038  SE(2,4r=-(CClTl.0)*A4+t7.0*B4+5.0*AB21*COEF 

00039  SKf  l.llssf  17.ft*frct  +  t  .0l*A4+fC9+t  ,  01*84*5. 0*AB2.1.*C.OEF .  . 

40 

41 

42 

00040  SK(3,4)s-( (C2+I .0)*B4+17.0#A4+5.0#AB2)*COEF 

00041  SK(4,4)a(l7.0*(Cl‘+1.0)*A44l7.0«(C2+1.0)*B4+25.0*AB2)*COEF 

O0042  iFfHTERM  _EQ-  41  GO  TO  20 

43 

44 

43 

00043  C 

00044  SK(2.5}a«*(16.O#(CC7Tl.0)»B4+A4'»'4,0-*AB2)*COEF 

00045  5K  (3*6  )=-f1  6  -0*  fCCl+1  .0}*A4*B4*4in *AB23*C0EF . 

46 

47 

48 

S0046  SK( 4#5)s(A4Tl6.0*B4+4.O*AB2)*COEF 

03047  SK(4,6)=(16.0*A4TB4+4,0*AB2)*COEF 

33048  SK(5.5)s(fCl+1.0)*A4+97.0#CCC2+1.0)*B4Tl1.0*AB2)*COEF 

49 

30 

31 

00049  SK{6»6)  =  C97.O»CCCt4i.0)#A4-»*(C2  +  1.0)*B44t3.0*AB2}*COEF 

03050  IFfHTERM  ,EQ.  6)  GO  To  20 

c  '  . : . . : . . 

52 

33 

54 

00352  SK(2»7)a(A4+l6.o*B4+4.0*A82)*COEF 

03053  SK(3,B)b(16.0*A4+b4+4v0*AB2)*COEF 

30054  SKf4.71s-M7.0*A4+16.0*(C2  +  1.0)*B4+20.0*AB2)#COEF 

53 

56 

3 

0O055  SK(4#8)=-(I6,0*(Ci+1.0)*A4+17.0*B4+20*0#AB2)*COEF 

30055  SK(5#7>*-(eCl+l*0)*A4  +  97.£j*644t3,0*AB2)*COEF 
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